Conversion of charge input into Kondo response 
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We show how the charge input signal applied to the gate electrode in double quantum dot may 
be converted to a pulse in the Kondo cotunneling current as a spin response of a nano-device in 
a strong Coulomb blockade regime. The stochastic component of the input signal comes as the 
infrared cutoff of Kondo transmission. 
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Current interest in charge-spin conversion effects is 
spurred by challenging prospects of spintronics. Most 
mechanisms of such a conversion are related to the in- 
terconnection between electrical and spin current due to 
spin-orbit interaction (SOI) [l|, which results in spin ac- 
cumulation near the sample edges. Such an accumulation 
in three- and two-dimensional electron gas in elemental 
and III-V semiconductors may result in spin-Hall effect 
Q and positive magnetoresistance Q . It was argued also 
that the Rashba-type SOI in a quantum dot induces spin 
current by a modulation of the voltages applied to the 
leads in a three-terminal device [4] . A spin Coulomb drag 
effect should be mentioned in this context, which results 
in spin polarization of the charge current due to intrinsic 
friction between electrons with different spin projections 
induced by Coulomb scattering 0, 0] . 

In all these propositions the possibilities of conversion 
of charge current into spin current were discussed. It is 
possible also to try to use the external electric field for 
the generation of spin current or another spin response. 
One such idea was formulated recently for light emitting 
diodes (LED) based on conjugated polymers 0, where 
dissociation of excitons in a strong enough electric field 
may result in the accumulation of up and down spin den- 
sities near the two ends of the LED. 

In this paper we show that the charge input signal 
applied to the gate electrode in a double quantum dot 
(DQD) may be converted into a pulse in the Kondo co- 
tunneling current, which is in fact the spin response of 
DQD under strong Coulomb blockade. We consider the 
mechanism of activation of internal spin degrees of free- 
dom by means of a time-dependent potential applied to 
an additional electrode to DQD in a contact with metal- 
lic leads in the so called T-shape geometry (Fig. [I]). The 
electrons in this device are described by the Hamiltonian 



H= E {ejnj + Qjn^+V^idUr. + H.c.ll) 

j—l,r;a <J 

+ Yl £ kb n kba + WY,i C kA° + h - C + Vg(t)nr 
b—s,d ka 

where the three first terms represent electrons in the 
DQD, electrons in the leads, and the lead-dot tunneling. 





FIG. 1: Left panel: DQD in a T-shape two-terminal geom- 
etry; right panel: evolution of energy levels in DQD as a 
function of scaling parameter £ = \nDo/D (see text). 



The last time-dependent term stands for the gate voltage 
applied to the dot r. Here = c\ ha c kh(J is the occu- 
pation number for band electrons with the wave vector 
k and spin a in the source and drain lead, n Jcr = dj a dj a 
the electron number in the left and right dot, and Qa r ) 
Coulomb blockade parameters. The tunneling Hamilto- 
nian involves only electrons in the left dot. Only the the 
even standing wave Ck a = (cksa + Ckda)/V2 enters the 
tunneling Hamiltonian. 

The time independent part H^ t of the dot Hamilto- 
nian may be diagonalized in the two-electron charge sec- 
tor M = 2, H$ t = J2a E a \A)(A\. The three lowest states 
at V g = are spin triplets \Tv) with spin projections 
v = 0,±1, a singlet 1*5) and a charge transfer singlet 
exciton \E) with energies Q 

E T = s l +s r ,Es = E T - 2(3V, E E = 2e l + Q l + 2/3V, (2) 

where f3 = V/Aes and £(i, r ) are the energy levels. The 
static part V g (0) of the gate voltage is incorporated in e r . 
Equations ([2]) are obtained for Q r ^> Qi, (3 <C 1. The 
ground state of an isolated DQD is a spin singlet Es- 

We study the influence of the gate voltage containing 
both coherent and stochastic components 



V g (t)-V g (0)=v g (t)=v g (t) + 5v g (t) 



(3) 



on the current through DQD. Here v g (t) = v g cos Fit is a 
coherent (deterministic) contribution and Sv g (t) a noise 



2 



component which is determined by its moments 



Sv g (t)=0, 5v g (t)5v g (t')=v2f(t-t') (4) 

The overline stands for the ensemble average and the 
characteristic function f(t — t') will be specified below. 
We incorporate V g (t) into the energy levels ([2|) by means 
of a canonical transformation 0, flo| 



H dot = t/^dotc/r 1 - ih -gf u l 



(5) 



with Ui = exp[— i<$>i(t)n r ] and the phase given by 

i r* 



*i(f) = \ J dt>v g (t>). 



(6) 



As a result, in the lowest orders in V g , the time dependent 
part of the dot Hamiltonian acquires the form 



5H dot (t) = -V ( + ) 



r(2) 



where 4 p) = Y.M**i* + (-1) p 4A*1. 

= § / / ^"/(*' - *")• 



(8) 



Thus, the time-dependent gate voltage induces coherent 
and stochastic interdot tunneling in DQD described by 
two terms in ([7j). One may formally include this term 
in the definition of the levels ([2]) and rewrite the dot 
Hamiltonian as Hdot(t) = X1a^a(^ AA where the level 
Et remains intact because the charge fluctuations do not 
influence spin degrees of freedom, whereas 



E s (t) = E s - S ad (t) - 5 atiS (t) 
E E (t) = E E + 5 ad (t)+5 stjE (t). 



(9) 



where the coherent part of renormalization 6 ad (t) = 
(2V 2 /AEs)®i(t) mav be considered adiabatically under 
the condition S^/Aes <C 1, whereas the stochastic cor- 
rection S s t y \(t) = (V 2 /4Aax) needs a more refined 
treatment. Here Aaa> = \E\ — E f A \. 

Coherent and stochastic components appear also in 
the cotunneling part of the effective Hamiltonian H cot 
which may be derived by means of the time-dependent 
Schrieffer- Wolff (SW) transformation |, [Io[. Unlike in 
the standard case of spin 1/2 quantum dots [sj, the SW 
transformation applied to DQD with the spectrum ([9]) 
intermixes the states |A). This intermixing is described 
by the operators |A)(A / |, which form together with diag- 
onal operators | A) (A | the set of generators of the SO (5) 
group. Ten generators are organized in three vectors S, 
P, M and one scalar A. Here S is the usual S=l spin 
operator, P and M are the vectors describing transitions 
between spin triplet \Tfi) and two singlets \S) and \E), 



respectively, A = —iS^/V2 intermixes the latter states 
under the constraint imposed by the Casimir operator 



(10) 



(see [10] for further details). The Hamiltonian H(t) — 
Hdot(t) +Hcot(t) m &y be rewritten in terms of the above 
group generators: 

H dot (t) = X - (^ T S 2 + ^P 2 +^M 2 ) -m(C-4) 
ffcotW = ^o T S-s + J 5 P-s + J E M-s. (11) 

The SO (5) symmetry is preserved by the last term in 
^dot by means of the Lagrange multiplier \i. Tilted cou- 
pling parameters are time-dependent. The charge-spin 
conversion mechanism under discussion is a manifesta- 
tion of the 50(5) dynamical symmetry of DQD. 

The Hamiltonian ([TT]) describes Kondo cotunneling 
through DQD in presence of time-dependent perturba- 
tions ([3|) Its coherent component is responsible for the 
conversion of the charge signal v g (t) into a Kondo- type 
zero bias anomaly in tunnel current response, which 
arises in the spin channel. Its stochastic component 
Sv g (t) results in the correction of charge noise into in- 
coherent corrections to Kondo cotunneling. 

The coherent part of the time-dependent Kon do p rob- 
lem can be solved in adiabatic approximation 0, [l(| . As 
a result of a time dependent Schrieffer- Wolff transfor- 
mation and elimination of high-energy states, additional 
renormalizations Ma(£) of the levels E\ arise where 
f = \n(D /D) is the scaling variable [1, [Tl| Then the 
singlet-triplet gap transforms into 



A st (i) = A° ST + M T -M S - Jad (*)• 



(12) 



and the singlet-triplet crossover (Fig. [TJ right panel) 
takes place, provided the self-consistent condition 
T K (A ST ,t) > A ST (t) is satisfied @, El, [H, [H, E3. 
Tk is a sharp function of Ast with a maximum Tko 
at Ast = (inset in the left panel of Fig. [2]). Its right 
slope is described by the ratio 



T K (A ST ,t) 



_A ST (t) 



(13) 



valid for intermediate asymptotic positive values of Ast 
at T K o/ A ST < 1 (dotted part of the curve T K (A)/T K (0) 
in the inset). Here rj < 1 is a universal constant. This 
sharp dependence is a key to the transformation mecha- 
nism of charge input into Kondo conductance response. 

We estimate the influence of v g (t) on the tunnel con- 
ductance G(T, t) at given T > Tk in a situation where 
the adiabatic temporal variations of Tk(Ast^) take 
place. Then the tunnel conductance obeys the law 



G/Go~ In" 2 (T/T K ). 



(14) 
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FIG. 2: (Color online) Left panel: Oscillations of ZBA in tun- 
nel conductance due to time- dependent gate voltage v g cos 
corresponding to the ST-gap oscillating around large positive 
A° ST > T K o (dashed line), around zero A down to T™ 71 = 
Tko/2 (dash-dotted line) and down to Tg in = (solid line), 
respectively. Right panel: time dependent Tk (notations are 
the same as for the left panel); insert: Tk as a function of 
Ast- Calculations are performed at T/Tko = 1.4, 77 = 0.5. 



Substituting ([T3j) in (jHj), one gets 

G(t)/G ~ (HT/T K0 ) - V ln(T K0 /A ST (t)r 2 (15) 

The curves G(Qt) shown in Fig. [2] (left panel) describe 
transformations of an oscillating signal in the charge per- 
turbation v g (t) into oscillations of Kondo-type zero bias 
anomaly (ZBA) of tunnel conductance. The transforma- 
tion effect is especially distinct provided the oscillations 
of v g (t) change the sign of Ast, i.e. induce an S — > T 
crossover (solid and dash-dotted curves in Fig. [2] corre- 
sponding to the solid part of Tr-(A) curve in the inset). 
It is worthwhile to notice that this mechanism of the 
adiabatic transformation of a charge signal into Kondo 
response is close to the one proposed for Kondo shut- 
tling [TBj where the source of time-dependence are the 
nanoelectromechanical oscillations of a quantum dot. 

The mechanism which converts a stochastic compo- 
nent Sv g (t) of the input signal into a stochastic spin re- 
sponse is quite unusual. Instead of dephasing due to 
time-dependent spin flip processes [9] , stochastization of 
the energy spectrum of DQD results in the loss of a Curie- 
type spin response at some characteristic energy (. This 
effect is related to the time dependence of the factor ji(t) 
in the Hamiltonian ([TT]) . Indeed, inserting ([9]) into ([TT]) . 
one may write the stochastic part of i^dot as 



rrst 
11 dot 



[5 st , s (t)P 2 - 5 st , E (t)M 2 ]/2 



(16) 



Unlike the adiabatic part of time dependent energy lev- 
els E\(t) incorporated in ([T2]h this term describes fluc- 
tuations due to the dynamical symmetry of DQD. At 
energies ~ Tk <C Ee the state \E) is frozen out. Then 
instead of the exact Casimir constraint ([TP]) for the SO (5) 
group, one deals with a fluctuating constraint C = S 2 +P 2 
for the reduced group 50(4) describing an ST multiplet 
where the fluctuating part may be written in the form 
u st (t)P 2 = -u st (t)S 2 with u st (t) = 6 atiS (t)/2. At T 
where the singlet is also frozen out, one arrives to the 



effective dot Hamiltonian 
1 



H dot (t) = -E T S 2 -fi ad (S 2 -2)-fi^ s (t)S 2 

*v=0,±l 



(17) 



where the fermion representation for S=l, \Tv)(Tv'\ = 
flf v t, is used in the second line. Here e = Er/2, and 
the time-dependent chemical potential for spin fermions 
is defined as fi(t) = /i a d — // s t,s(t). The stochastic com- 
ponent of u may be treated as a random potential in 
the time domain which describes the fluctuations of the 
global fermionic constraint [l0[. Then the propagation 
of spin fermions in the random time-dependent potential 
/i s t,s(£) m &y be studied by means of the "cross technique" 
developed for the calculation of electron propagation in 
a field of impurities randomly distributed in real space. 

The frequency £1 in v g (t) is the slowest frequency in 
our problem, and /x s t,s(t) also varies slowly in time, so 
that the relaxation time r = h/j in the noise correlation 
function D(t - if) = h 2 ~ exp[~7(t - t')\ is a 
longest time in the model. Then one may take the limit 

2C 2 7 



D(lu) = lim 

7^0 UJ 



2 +7 2 



(18) 



for its Fourier transform, £ being the bandwidth of the 
Gaussian correlation. In this limit the averaged spin 
propagator describes the ensemble of states with chem- 
ical potential /i s t= const in a given state, but this con- 
stant is random in each realization [l6[. Thus the prob- 
lem of decoherence of the spin state in stochastically per- 
turbed DQD is mapped on the so-called Keldysh model 
I?! 0, [l9[ originally formulated for (5-correlated impu- 
rity scattering potentials in momentum space. The prob- 
lem can be solved exactly and the decoherence time is 
fixed by £ of the Gaussian. 
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FIG. 3: Feynmann diagrams for the self energy E with ver- 
tex corrections, triangular vertex T, spin susceptibility \ and 
parquet Kondo loops K (from left to right). Solid, dashed 
and wavy lines denote bare spin- fermion propagator G, con- 
duction electron propagator g and correlation function D. 

The solution of the Dyson equation obeyed by 
G~ 1 (e) = e + fi ~ S(e) (Fig. 0) [20], the Fourier 
transform of the spin- fermion propagator G^ u (t — t') = 
UMtt(f))n = -i([/„(t)/J(*')] + >, is given by 



m - / 



duo 
2^ 



r(e, e - u\ uo)G(e - uo)D(uo) = C 2 r(e, e; O)G(e) 



(the index v is omitted, since the fluctuations of u 
are related to the global U(l) symmetry). Using the 



Ward identity for a triangular vertex shown in Fig. O 
r(e, e;0) = dG~ 1 (e)/de, we transform the Dyson equa- 
tion into differential equation 



( 2 dG/dx + xG 
where x = e + /ig with a solution [l 
1 



1 = 



G K (x) = 



CV27T 



id 



(19) 



(20) 



representing the set of spin states with stochastic chemi- 
cal potential averaged with a gaussian exponent charac- 
terized by the dispersion (. Remarkably, G R (x) has no 
poles, singularities or branch cuts. 

To check the spin properties of stochasticized DQD, 
we calculate its spin response at finite T determined as 
x(iu m ) = T ^ n QiiuJm + ie n )Q{ie n )T(ie n , ie n + i(x} m \ iu m ) 
where Q(ie n ) is the Matsubara continuation of ([20]) on to 
the imaginary axis (see Fig. [3|). The same Ward identity 
provides the exact equation for the vertex 



iR\2 



r w (e,e;0) = (eG K -l)/C^(G*) 
which allows to calculate the static susceptibility 
1 



X(" = 0,T) 



8ttC 



dxe x / 2 £tanh 



2T 



(21) 



(22) 



From its temperature dependence plotted in Fig. (4) we 
see that the the Curie- type spin response x(0, T) ~ 1/T 
at high T ^> ( transforms into a constant at zero T, 
x(0,0) ~ 1/C- This means that the DQD looses at 
{cj, T} <C C the characteristics of a localized spin due 
stochastization, so it cannot serve as a source of Kondo 
screening at low energies. The frequency dependence 
of Imx(c<;, T) (Fig. 4, inset) confirms this conclusion. 
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FIG. 4: (Color online) Static susceptibility x(0,T). Inset 
shows a frequency dependence of Imx with a maximum at 

Such a behavior of the DQD affects the Kondo response 
described by H cot ([TT]) . Time dependence in the ver- 
tex T (Fig . [3]) results in dephasing of Kondo cotun- 
neling [lOj], which, however, is exponentially weak at 
T <C A St- More significant is the influence of stochas- 
tization of the DQD. Inserting ([20]) into the Kondo loop 
K(ie n ) ~ J 2 T^2 m J j0^g(p,iuj rn )G(ie n ±iu m ) respon- 
sible for logarithmic singularity in conventional Kondo 
scattering (Fig. [3]), one obtains a combination of log- 
arithmic, hypergeometric and imaginary error functions 



:PoJ 



1*1 



3 o H 

2' ' 2C 2 



7rErf i 



V2C, 



where 7 = InC is the Euler constant. In two limiting 
cases of low and high temperatures relative to the dis- 
persion ( of the noise spectrum, it leads to following ex- 
pressions 

( poJHD/T), T»C 
K(e -> 0)/ J = { (23) 
{ p Q J\n{D/Q, C»r 

Thus the noise amplitude plays the role of the infrared 
cut-off (similarly to Kondo-spin glass problem [ll|). This 
cut-off distorts the coherent pulses in ZBA (Fig. [2]) pro- 
vided is comparable with Tk • 

In conclusion, we presented and discussed the con- 
version mechanism of a time-dependent coherent and 
stochastic input signal at a gate electrode into a Kondo- 
type spin response in tunnel conductance. This mecha- 
nism is related to the dynamical symmetry of DQD. 

Authors are indebted to F.Marquardt, N.Prokof'ev, 
M.V. Sadovskii and T. Ziman for valuable advices. 
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